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Abstract 

We consider a model of a distribution feeder connect- 
ing multiple loads to the sub-station. Voltage is con- 
trolled directly at the head of the line ( sub-station), how- 
ever, voltage anywhere further down the line is subject 
to fluctuations, caused by irregularities of real and re- 
active distributed power consumption/generation. The 
lack of a direct control of voltage along the line may 
result in the voltage instability, also called voltage col- 
lapse - phenomenon well known and documented in the 
power engineering literature. Motivated by emerging 
photo-voltaic technology, which brings a new source 
of renewable generation but also contributes significant 
increase in power fiow fiuctuations, we reexamine the 
phenomenon of voltage stability and collapse. In the 
limit where the number of consumers is large and spa- 
tial variations in power flows are smooth functions of 
position along the feeder, we derive a set of the power 
flow Ordinary Differential Equations ( ODE), verify phe- 
nomenon of voltage collapse, and study the effect of dis- 
order and irregularity in injection and consumption on 
the voltage proflle by simulating the stochastic ODE. We 
observe that disorder leads to nonlinear ampliflcation of 
the voltage variations at the end of the line as the point 
of voltage collapse is approached. We also flnd that the 
disorder, when correlated on a scale sufficiently small 
compared to the length of the line, self-averages, i.e. the 
voltage profile remains spatially smooth for any individ- 
ual realization of the disorder and is correlated only at 
scales comparable to the length of the line. Finally, we 
explain why the integrated effect of disorder on the volt- 
age at the end of the line cannot be described within a 
naive one-generator-one-load model. 



1 Introduction 

Voltage collapse/instability is a severe disturbance but 
also one of the most interesting and most discussed non- 
linear phenomenon in power engineering, with the bibli- 
ography [ 1 1 published in 1998 containing 308 references 
and a review in 2000 |2 | containing 132 references. Ac- 
cording to 1 IJ, the first paper related to voltage instabil- 
ity appeared in O, and the first criteria for detecting the 
point of voltage collapse was proposed in [Tl. A number 
of comprehensive books |51 |6l 13 and a book chapter in 
lIHl are written on the subject. 

In this manuscript, we give an applied math and sta- 
tistical physics prospective on description of the voltage 
collapse phenomenon in the spatially continuous (ODE) 
limit, and we also begin to investigate effects of spa- 
tial disorder in power injection/generation on the volt- 
age instability. To the best of our knowledge, an ODE 
description of power flows was discussed in the power- 
engineering literature only in the context of the electro- 
mechanical waves 191 do) 1 and not to describe voltage 
collapse and instability. We are also not aware of any 
discussion in the literature of the effect of spatial disor- 
der, and resulting structural stability, in power consump- 
tion along a power line[^ 

A natural question that should be answered before 
the development of ODE models is whether the contin- 
uous models are useful for modeling discrete systems 
like the power grid. Similar questions have been an- 
swered multiple times in the field of statistical mechan- 
ics which aims at establishing connections between mi- 

' Note, that the effect of spatio-temporal noise in consumption on 
the voltage collapse in a transmission system was already discussed in 
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croscopic models of various materials (for instance gas 
of interacting particles) and their thermodynamic prop- 
erties. The main lesson learned from decades of research 
in this area is that all details of microstructure are in fact 
not that important when one is focused on describing 
the large-scale properties of the system. The starting 
point of our approach is the observation that power net- 
works share many properties with statistical mechanics 
systems. Huge number of individual loads contribute 
to the power flow; however, modification, inclusion or 
removal of any individual load will not lead to any vis- 
ible macroscopic change. This general macroscopic ap- 
proach motivates us to ask the question: is a coarse- 
grained, continuous ODE description a proper tool to 
describe macroscopic phenomena in power grids? 

To answer this fundamental question we extended the 
ODE analysis of the voltage collapse/instability to an 
inhomogeneous/disordered case. Within the ODE de- 
scription, we consider effects of temporally frozen but 
spatially varying irregularity and disorder in real and re- 
active consumption and generation on the voltage pro- 
file in a representative feeder line of a power distribution 
system, rj Assuming that the correlation scale of the dis- 
order is smaller than the length of the feeder (but larger 
than the inter-distance between neighboring nodes), we 
study the resulting stochastic power plow (PF) ODEs in 
simulations. Our main findings consists in the following 
three observations: 

• The effect of disorder on the voltage profile, and 
specifically in voltage variations at the end of the 
line, is amplified in a significant way with the pa- 
rameters of the base solution (uniform component 
of consumption along the line) as the point of the 
voltage collapse is approached. The probability of 
observing voltage collapse is increased in a nonlin- 
ear fashion in the presence of disorder. 

• Any individual configuration of the disorder results 
in a smooth spatial profile of voltage and power 
flows, correlated at the scale roughly associated 
with the length of the feeder line, and which is cer- 
tainly much longer then the correlation scale of the 
disorder We attribute this self-averaging effect to 
the fact that solutions for the power flows, but es- 
pecially for the voltage, emerge as spatial integrals. 
The observed smoothness of the solution implic- 
itly justifies using the reduced ODE approach in the 
power systems. 

• The effect of disorder on the voltage at the end 
of the long feeder line cannot be described prop- 

-Note, however, that the developed formahsm is generic and as 
such it as well applies to linear segments of transmission lines, which 
will also be discussed biiefly in the following. 



erly by a reduced one-generator-one-load model, 
i.e. the effects revealed by the stochastic ODE de- 
scription cannot realized in a lumped model by in- 
tegrating the injection and consumption of real and 
reactive powers along the line. 

The material in the manuscript is organized as fol- 
lows. Section|2]contains a technical introduction resolv- 
ing in the ODE formulation of power flows for linear 
segments of the grid. Section |3] analyzes phenomenon 
of voltage collapse and instability within the ODE ap- 
proach, discussing both a linear element of a transmis- 
sion system and a linear feeder line of a distribution sys- 
tem. Section|4]is devoted of analysis of a structural sta- 
bility of a feeder line subjected to static variations in 
consumption and injection along the line which are cor- 
related on a short length scale. Section|5]is reserved for 
conclusions and brief discussions of the path forward. 

2 Power Flows in Discrete and 
Continuous (ODE) Formulations 

Static Power Flow (PF) equations over a power 
graph/network, Q — (1^,2;), are reformulations of the 
AC Kirchoff laws for currents and voltages defined at 
any vertex of flie graph, a e 1^: 

= /5„ - 2^ VaVb sin(9„ - 0;,) . . 

-i:(v^v„v,cos(e„-e.))^^^' (1) 

i~a '^ab^^ab 

^« - E (va-VflViCos(e„ - e,,)) 2^'^ 2 

br^a '^ab^^ab 
b^a '^ab^^ah 

where /5„ , , v'„ , 0„ are real and reactive power injections 
(negative for consumption), potential, and phase at the 
node a E 1^, respectively. Here, b ^ a indicates that b 
and a are graph-neighbors connected by an edge {a,b) 
of E characterized by resistance R^h and inductance Xah- 
In transmission (high voltage) lines, resistance is signif- 
icantly smaller than reactance (Rab <C Xah), while in the 
distribution (low voltage) networks resistance and reac- 
tance are usually of the same order. 

Relations between the electric potential (voltage and 
phase) and powers (real and reactive) are universal, how- 
ever, what constitutes input or output vary for nodes of 
different types. Standard generator nodes are kept at 
constant voltage (assumed unity in the rescaled units), 
and thus their input pair is real power and voltage, i.e. 
they are {p,v) nodes. Regular consumer nodes draw 
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prescribed powers, therefore these are {p, q) nodes with 
real and reactive powers being the two input parameters. 
The largest generator of a system, or a node connect- 
ing the system to a bigger power network, is taken as 
a phase reference where both voltage and phase are as- 
sumed fixed, i.e. a (v, 9) node. Summarizing, the power 
flow equations Eq. \\\2\ constitute a system of algebraic 
equations expressing unknown (output) variables, which 
are real and reactive powers at the super node, reactive 
powers and phases at the generator nodes, and phases 
and voltages at the consumer nodes, via the set of known 
(input) node variables. 



Note that the {p,q) model with constant and time- 
steady p and <7 is a relatively crude approximation. In 
reality, p and q do depend on instantaneous voltage at 
the load v, the local rate of phase variation 9 (i.e. on 
the frequency), and on the rate of voltage variation v. 
Modeling the so-called static and dynamic 

[(/?,q')(9, v)] loads are classic and also difficult subjects 
in power engineering. See a (relatively) early discus- 
sions of the load modeling in HIl [131 |l4J and [15,, J6J 
for description of the current status quo and problems in 
the reduced but (still multi-parametric) composite mod- 
eling of loads within Western Electricity Coordinating 
Council (WECC) ^ presented as a part of the GE Pos- 
itive Sequence Load Flow (PSLF) software fl8|. The 
main difficulty in the microscopic modeling of loads 
comes from the uncertainty and variability of loads over 
time. Indeed, an aggregated load representing a node of 
the grid may consist of a very special combination of 
nonlinear and dynamic elements(like rotors) and non- 
linear and linear static elements (like resistors or induc- 
tances). On the time scale of interest here (seconds or 
minutes), an aggregated load load model is not expected 
to change in time, i.e. it is quenched. The situation 
poses the question of "learning" the nonlinear and dy- 
namic (but reduced) models of loads. 



It is common in a transmission network to have a 
linear segment connecting two generators with multiple 
loads drawing power from the line in between. Feeder 
lines are the most typical elements of distribution net- 
works where strong connection to the transmission net- 
work resides at the head of the line, and the line ends at 
the last load of the linear segment. Assuming that the 
loads are distributed uniformly and frequently along the 
line so that the inter-load spacing a is sufficiently small 
{a ^ L), the voltages and phases can be expanded in 
a Taylor series in a, v„+i = v„ + (3,v')„fl + (3^v)„a^/2 + 
0(fl^),9„+i =9„ + (a,9)„fl + (3?9j„flV2 + C(a3). Sub- 
stituting the expansions into Eqs. ( |l|2[ i and keeping only 



the leading 0{a^) terms, one arrives at 



= p,nech + {v%Q) +gv (a?v - V (3,9)2) , (3) 
= q„,ech + Pv (dh - V (3,9)2) - gd, {v%Q) , (4) 



where pmech = 0/a,qmech = q/a are real and reactive 
power densities]^ g = aR/{R^+X^) and ^ = aX/{R^ + 
X^) are conductance and susceptance measured in the 
units of the inter-load spacing, a. Note that the two 
second-order ODEs P|4| i can be restated as the follow- 
ing four first-order ODEs (written at p = 1 for simplicity 
of notations ) 



V w ; 



s 

w 

p+m 



\ 



(5) 



V ..3 (l+,?2)v / 



Flows of real and reactive powers 

P = -pv23,9 -gv3,v, Q = -pv3,v + gv23-9, (6) 

are generally dependent on the position along the line z. 
These two objects, which appears on the rhs of Eqs. ([3]l, 
are important characteristics of the flows. (See Ap- 
pendix|A]for more details.) 

For a transmission line, the four natural boundary 
conditions (for two second order differential equations) 
are known voltage at the two ends of the line [v(0) = 
v(L) = 1], zero (reference) phase at the head of the line 
[9(0) = 0], and a fixed value of the real power flow at 
the head of the line [/"(O)]. For a distribution line, the 
four natural boundary conditions become: known volt- 
age at the head of the line [v(0) — 1], zero (reference) 
phase at the head of the line [9(0) = 0], and zero flow 
(flux) for both real and reactive powers at the end of the 
line [P(L) = e(L)=0]. 

Note that the PF equations, in discrete (T]2i or con- 
tinuous ( |3|4[ ) forms, may show multiple solutions or 
have no solutions. It is also important to mention, even 
though we do not discuss in the manuscript, that the 
ODE-based static power flow framework naturally fits 
in a more general dynamic framework, accounting for 
dynamics and control of generators and loads. 



The subscript "mech" in p and q indicates tliat the powers are 
mechanical powers at the buses. In the following we will skip the 
subscript to lighten notations. 
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Figure 1: Lowest voltage along the [0;L] segment vs the 
length of the segment for fixed values of load consump- 
tion (p ~ —5 and ^ = — 1), conductance (g = 0.5), and 
susceptance (p = 1). The value of the real power in- 
jected in the line at L = is varied. Dashed curves show 
the respective voltage profiles along the critical segment, 
[0;L,]. 
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Figure 2: Profile of real (solid) and reactive (dashed) 
powers flowing along the line. The curves correspond 
to the critical segments with the parameters (and color- 
coding) from Fig. ([TJ. 



Profile at the point of Voltage Collapse 




position along the line, Z6[0;L] 



3 Static Analysis of the Continuous 
Model in the Spatially Uniform 

Case 

Our first step is to verify that transition to the ODE de- 
scription of a linear segment (for a transmission network 
or a feeder line) preserves the phenomenon of the volt- 
age collapse/instability observed in discrete power sys- 
tems. 

We consider the ODE PF model ( |3|4| i of a transmis- 
sion line with constant p, q, g, and p homogeneously 
distributed along the line z G [0;L] with voltage fixed at 
the both ends of the line v(0) = v{L) = 1. We also fix 
the value of a real power P{z = 0) injected at the head of 
the line, and we study the dependence of the voltage and 
phase profile along the line as a function of L, V9(0), p 
and q. The set of two ODEs allows analytic solution (in 
quadratures) only in the special (but unphysical) case of 



Figure 3: Dependence of the voltage on the density of 
real power consumption along the feeder for L = 0.5 and 
3,-9(L) = 3j:V'(L) = 0. Doted and dashed curves corre- 
spond to stable and unstable branches respectively. 



q=p/5, g=p-1 , v(0)=1 , P(L)-Q(L)=0 




Density of Consumed Real Power, -p 



p = 0, therefore we largely rely here on a straightfor- 
ward numerical study of Eqs. (3]4i. Simulations are 
performed using the shooting method. 

Typical results of direct illustrating phenomenon of 
voltage collapse are shown in Figs. (T]2i. In this static 
interpretation, voltage collapse consists in the disappear- 
ance of a valid solution for sufficiently long lines of 
length L > Lc, where thus L^. is the critical value cor- 
responding to the point of voltage collapse. For L < L^, 
valid solutions appear in pairs, one stable and one un- 
stable correspondoing to high and low voltage, respec- 
tively. For L < Lc, these two solutions vary with L, and 
they merge at the critical point L — L^. Modifying g/^ 
does not change the qualitative shape of the curves in the 
Figs. (T]2i. We note a non-monotonic dependence of Lc 
and Vc (voltage at the point of collapse) on P(0). 

Main conclusion one draws from Figs. ( |l|2| i and re- 
lated analysis is that the "nose shape" behavior is a uni- 
versal phenomenon. The stable solution corresponds 
to higher throughput (of both real and reactive) power 
over the line. On the quantitative side, position of the 
nose (both in terms of the critical length of the line and 
the critical voltage) is a non monotonic function of the 
line parameters: the power consumed along the line, the 
power dissipated within the line, and power transferred 
through the line. The line is the longest in the symmetric 
situation of the zero throughput. 

Next, we analyze a feeder line of fixed length L where 
the end of the line is hanging free and the amount of 
real and reactive powers injected at the beginning of the 
feeder is adjusted to the following four boundary condi- 
tions, v(0) = 1,9(0) = 0,P(L) = Q{L) = 0. The bound- 
ary value problem is studied with the help of the bvp4c 
Matlab solver. The results are shown in the Figs. ( 3]4 i 
and describe the "nose shape" for a feeder of the fixed 
length. 

Here we conclude that the two solutions, one stable 
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Figure 4: Dependence of the real and reactive power at 
the head of the Hne on the density of real power con- 
sumption, complementing description of Fig. Q. Solid 
and dashed curves show stable and unstable branches re- 
spectively. 



q.pffi, a.p.1,v(0|.1, P(L).Q(L|.0 




Figure 6: Multiple (twenty trials per case/color) configu- 
rations of the disorder resulted in voltage profiles shown 
in Fig. (j5]). 
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Figure 5: Voltage profile (twenty trials per case/color) 
for three base configurations of the feeder line and added 
Gaussian finite-correlated disorder (in both p and q) of 
the same strength. (Red blue and green corresponds to 
p = —2.5,-3 and —3.5 respectively.) Low and Upper 
branches (solid and dashed) correspond to stable and un- 
stable solutions. See text for further details. 
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and another unstable, are observed for a feeder of a fixed 
length when the amount of real and reactive power con- 
sumed (homogeneously) through the feeder is smaller 
than the critical value, corresponding to the nose of the 
curve in Fig. Q. We also observe the power injected 
at the head of the line is smaller for the stable solution 
than for the unstable solution, confirming the intuition 
that the choice of the stable solution is consistent with 
the minimal energy/power principe. 



4 Structural Stability Analysis 

Structural stability pertains to modification of the typi- 
cal"nose" curve shown in Fig. ([3]l because of structural 
disorder in the p and q, i.e. inhomogeneity in the p,q 
distribution along the line. 

One important phenomenon, illustrated in Fig. (j5]l. 



Figure 7: Comparison of the actual disorder effect 
(one realization from the ensemble shown in Fig. (|6]l) 
on the voltage profile (star-line) with an aggregated 
effect (dash-line), where in the latter case p and q 
along the line are kept uniform and correspondent to 
PQ and qo with added small values, /q dp{z')dz' /L and 
jQdg{z!)dz' /L respectively. The green is not seen be- 
hind red - the curves virtually coincide. 

Effect of Structural Disorder on feeder line: L=0.5, o=0.5, t=0.01 , q=p/5, p=g=1 , p=-3.5 
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Figure 8: Configuration of disorder correspondent to the 
profiles shown in Fig. (j7]i. 



An instance of ctisortier: c=0.5, t^O.OI , N=1 000 
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is the amplification of the effect of disorder when ap- 
proaching the point of vohage collapse. We tested 
three sample cases from Figs. P|4[ ) corresponding to 
an average density of consumption po = —2.5; —3; —3.5 
with disorder in p,q added. The disorder, 5p(z) and 
5^(z) where p{z) = PQ + 5p(z) and q{z) =qo + 5q{z), 
is modeled as zero-mean, colored Gaussian noise char- 
acterized by (5p(z)5p(z')) = 2a^exp(— |z — z'l/x) and 

(5^(z)5,(z')> - 2(o/5)2exp(-|z-z'|A). where o = 
0.5, T — 0.01. Multiple realizations of the disorder are 
shown superimposed in Fig. (j6]). 

As illustrated in Figs. ( |7|8[ l, the shift in the voltage 
profile due to disorder is not associated with a sim- 
ple renormalization of the cumulative loads / p{z!)dz! /L 
and / q{{)d{ jL along the line. The fact that the curves 
without disorder and uniformly distributed disorder co- 
incide in Fig. (|7]l is incidental. The order of curves and 
relative difference between them may vary from one dis- 
order configuration to another, but in general differences 
between three type of curves (no-disorder, disorder, and 
uniformly distributed disorder) are of the same order. In 
some cases the disorder may significantly degrade the 
situation but in others it may improve it. In some rare 
cases the correction added by disorder was super-critical 
thus resulting in the algorithm failing to find any solu- 
tion. 

Two main conclusions drawn from these simulations 
are as follows. 

• In spite of the the fact that the amount of disorder 
in p and q added was identical in all the three cases, 
the spread of voltage was nonlinearly amplified as 
criticality is approached (the point of voltage col- 
lapse). The nonlinear amplification is associated 
with the fact that the closer the system is to the 
point of the voltage collapse, the easier it becomes 
for disorder to push the system even closer towards 
instability or possibly overshoot into the unstable 
regime. Since the behavior is highly nonlinear in 
the vicinity of the point of collapse, the effect of 
amplification is nonlinear too. 

• The effect of disorder on the voltage profile self- 
averages: the voltage profile decreases monoton- 
ically even after the sub-critical configuration of 
the disorder is added (when the system still has a 
dynamically stable solution). Moreover, we have 
also observed that the overall effect of the disor- 
der cannot be explained simply as a renormaliza- 
tion of the total power consumed along the line. 
Adding disorder may push the system to criticality 
or relieve it somewhat. The self-averaging feature 
of the resulting voltage profile is associated with 
the nature of power flow equations and boundary 



conditions. Indeed, the DistFlow ODE representa- 
tion (explained in the Appendix) makes this clear. 
Eqs. Q are integrated backwards in z from L, sub- 
stituted into Eq. ( lOl, and integrated forward in z. 
This double-integral structure results in smoothing 
of the input/disorder in the voltage profile/output. 
It is also obvious that the resulting profile is not a 
simple integral of the disorder, but rather a nontriv- 
ial convolution of the disorder with a kernel that is 
nonuniform along the feeder line. 

5 Conclusions and Path Forward 

In this manuscript we developed an ODE approach 
to static analysis of voltage instability collapse in lin- 
ear segments of transmission and distribution power 
systems. The technique allows the use of powerful 
ODEs/PDEs tools of applied mathematics and statistical 
physics for deriving a reduced, coarse-grained descrip- 
tion of power flows and voltage variations along these 
lines. The technique was validated on examples of lines 
with uniform and homogenenous distributions of con- 
sumption/generation of real and reactive power. The PF 
ODEs shows voltage collapse and instability, similar to 
these observed in finite systems. Then we used the tech- 
nique to study effects of variations and disorder in con- 
sumption/generation on voltage profile. We observe that 
the effect of disorder, in particular on the value of volt- 
age at the end of the feeder line, is amplified in proxim- 
ity of the voltage collapse. Moreover, we observe that 
the resulting voltage profile remains smooth even un- 
der action of a short-correlated disorder, which justifies 
the use of the reduced ODE approach for power system 
modeling. 

This manuscript is our first paper reporting results of 
the ODE analysis of voltage stability in power systems. 
We envision extending this research along the following 
(and possibly many other) directions. 

• Introduction of dynamic (deterministic and 
stochastic) modeling of load and generation and 
studying dynamic phenomena underlying the static 
analysis of this paper. We will be developing 
practical tools for analyzing regions of dynamic 
stability in linear segments of power systems 
within the reduced ODE/PDE approach. We 
anticipate using this approach to extend the prob- 
abilistic distance to failure methods of lfT9l l20ll 
based on large deviation techniques to account 
for dynamic failures associated with the voltage 
collapse/instability. 

• Our ODE approach can be extended to analyze 
two-dimensional regions of power grids, such as 
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coarse-grained models of densely populated areas 
of the Eastern Interconnection in US. 

• The newly developed stochastic ODE/PDE tech- 
nique will be an instrument for developing new 
control schemes based on detection of dangerous 
stochastic contingencies (see aforementioned dis- 
cussion of the distance to failure approach) be- 
fore they actually happen and mitigating them ef- 
ficiently by using a variety of available energy and 
power resources. 

• We plan to use this ODE/PDE approach to study 
other types of instabilities in power grids, in 
particular these associated with synchronization- 
desynchronization transitions in phase dynamics 
ETl . effects of inverse cascade (broadening of 
the desynchronization frequency spectra) [22 J . and 
electro-mechanical waves ||9l- 

• It will also be important to extend the ODE/PDE 
model-reduction ideas to extracting/learning static 
and dynamic model parameters (especially these 
related to small scale distributed generation and 
loads) from measurements, in particular com- 
ing from the emergent Phasor Measurement Unit 
(PMU) technology L23J. 
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A DistFlow Representation 

In the case of a linear segment it may be convenient to 
restate the power flow equations in terms of the power 
flowing at any point of the segment. In accordance with 
the original (spatially discrete) formulation of Eq. ( |l|2| i, 
one establishes the following relation between the power 
leaving a site n and moving towards the next site along 
the line, « + 1, and respective characteristics of the po- 



tential: 



X 



R 



v„v„+isin(e„-e„+i) 



Q" = ^^2~|:^^n(^n-^«+icos(e„-e„+i)) 

In the continuous limit these relations become Eqs. (|6|, 
while Eqs. (3]4i turn into 



pZ iQZ pZ, QZ 

= p-d,P-r 0^q-d,Q-x (9) 

where r — R/a, x — X /a are resistance and inductance 
densities, and, in accordance with Eq. (j6]), P and Q are 
related to v as follows 



rP+xQ 



(10) 



In their discrete version the set of Eqs. ( |9|10| l for a line 
are called DistFlow equations ll24l |25]| . (See also 
.) 
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